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IN***' 


I.  INTRODUCTION 

Most  models  for  simulating  the  concentration  of  obscurant  in  a  smoke  or 

dust  cloud  make  use  of  gaussian  plumes. ^  ^  This  concentration  must  be  in¬ 
tegrated  along  some  line-of-sight  that  connects  an  observer  to  a  target.  The 
complexity  of  this  integration  usually  dictates  that  it  be  performed  numeri¬ 
cally  at  great  cost  in  terms  of  both  computer  memory  and  execution  time.  In 
this  report,  an  approximate  method  for  calculating  this  concentrat ion- 
pathlength  in  closed  form  is  described. 

This  investigation  is  part  of  an  ongoing  examination  of  the  effect  of 
battlefield  obscuration  on  the  performance  of  missile  (and  selected  nonmissile) 
v/eapon  systems  as  part  of  the  Concepts  Analysis  and  Validation  work  area  of 
the  A214  Missile  Technology  Program.  The  results  of  this  examination  will  be 
used  in  the  formulation,  analysis,  and  evaluation  of  present  and  conceptual 
missile  weapon  systems. 


1 1 .  CONCENTRATION  -PATHLENGTH  DEVELOPMENT 

The  transmission  through  an  obscurant  from  an  observer  to  a  target  (point 
in  space  to  point  in  space)  is  commonly  expressed  as 

T  =  exp  [-aC^ ]  (1) 

where:  T  =  transmission, 

2 

a  -  extinction  coefficient  of  the  obscurant  (cm  /g) ,  and 

2 

C^  =  concentration-pathlength  through  the  obscurant  (g/cm  ). 


The  concentration-pathlength  is  usually  given  by 
target 

C^=  /  C (x,y , z)  dt(x,y,z)  (2) 

observer 

where:  C(x,y,z)  =  concentration  of  obscurant  at  point  with  coordinates 

3 

x,y,z  (g/cm  ),  and 

£(x,y,z)  =  line-of-sight  (cm). 

The  most  common  coordinate  system  used  is  a  cartesian  coordinate  system 
with  origin  at  the  formation  point  of  the  obscurant  cloud  (the  cloud  co¬ 
ordinate  system).  The  x  axis  of  this  system  lies  in  the  ground  plane  along 
the  wind  direction,  positive  downwind.  The  z  axis  is  vertical  and  the  y  axis 
is  oriented  so  that  the  system  is  right-handed.  An  alternate  coordinate 
system,  called  the  cloud  centroid  coordinate  system,  has  the  same  orienta¬ 
tion,  but  the  origin  is  at  the  cloud  centroid.  This  coordinate  system  is 
time  dependent. 


1 


The  observer  and  target  have  coordinates  x  ,  y  ,  z  and  xfc ,  \r  ,  z  in  the 

o  o  o  t  t  t 

cloud  coordinate  system.  The  centroid  of  the  cloud  has  coordinates  uCt-t^), 

0,  z(t-tQ),  where  u  is  the  ground  plane  wind  speed,  t^is  the  initial  formation 
time  of  the  obscurant,  and  z(t)  is  the  centroid  rise  function. 

The  concentration  of  obscurant  is  most  commonly  expressed  as  some  form  of 
trivariate  gaussian.  One  of  the  simplest  of  these  which  contains  the  com¬ 
plexity  needed  for  this  report  is  that  of  the  Smoke  Effectiveness  Manual 
Model  (SEMM).^O  This  model  is 
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where : 

Q  =  munition  load  of  ^ cur  ant  producing  material  (g)  , 

X  =  efficiency  factor, 
ft  =  yield  factor, 

a  ,  a  ,  a  =  cloud  standard  deviations  (time  dependent)  (m) , 
x  y  z 

x,  y,  z  =  cartesian  coordinates  of  points  in  cloud  coordinate  system  (m) , 
u  =  wind  speed  (m/s),  and 
t  =  time  since  cloud  formation  (s). 

The  cloud  rise  function  is  zero.  Equation  3  is  appropriate  for  an  instanta¬ 
neous  cloud. 

Two  methods  are  commonly  used  to  develop  the  explicit  form  of  the 
concentration-pathlength:  line-of-sight  (LOS)  parametr izat ion  and  coordinate 

system  rotation.  Each  shall  be  reviewed  separately.  In  LOS  paramet rization , 
the  coordinates  along  the  LOS  are  written  as 


x  =  x  -l-  Axn 
o 

(A) 

y  =  yQ  +  Ayn 

(5) 

Z  =  Z  +  Azn 
o 

(6) 

where : 

Ax  _ 

d 

(7) 

( v  -  v  X 

A  v  =  t  -  o 

c 

(8) 

(z  -  2  ) 

_  t  O 

(9) 

Tin*  distance  between  the  observer  and  the  target  is 


2 


d  =  [(xt'Xo)2  +  (Vyo)2  +  (zt-Zo)2] 


2  I  1/2 


(10) 


The  parameter  n  takes  on  all  values  between  0  and  d.  Equations  4  through  6 
are  substituted  into  Equation  3  to  yield 
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after  some  minor  rearrangement.  The  explicit  form  of  the  concentration- 
pathlength  is  then 


ci  = 


d 

/ 


C(n)dn 


(12) 


Although  only  the  last  two  terms  of  the  right  hand  side  of  Equation  11 
are  n  dependent,  the  integral  of  Equation  12  must  either  be  performed 
numerically,  the  square  of  the  exponent  completed  and  the  result  extracted 
from  a  table,  or  if  d  is  sufficiently  large,  an  approximation  that  d  =  «>  may 
be  introduced  and  the  integral  performed  analytically . H 


The  technique  of  coordinate  system  rotation  first  translates  into  cloud 
centroid  coordinates,  and  then  defines  the  scaled  coordinate  system  given 
by 


a 
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y' 


=  JL- 

a 
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(13) 


where  x,  y,  z  are  in  cloud  centroid  coordinates.  Thus,  the  observer  co¬ 
ordinates  in  this  system  are 
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(x  -  ut) 

t  O 


X 

o  o 


yo  =  T 


(14) 


o  a 


and  the  concentration,  Equation  3,  is 


C(x  ,y  ,z') 


2QXfi 


(2tt)3,/2o  a  a 
x  y  z 


exp  [-l/2jx’2  +  y’2  +  z'2j] 


(15) 


Next,  the  ground  plane  distance  (scaled)  between  observer  and  target  is 
defined  as 


d-  =[(x;-x;>2  +  (y^-y0)2] 

and  two  angles 


1/2 


*  arc  tangent 


and 


0  -  arc  tangent 


(yt_yo) 


df 


These  two  angles  define  a  rotation  matrix 

cos  (0)  cos  (cj>)  cos  (0)  sin  (4>) 

R  =  -sin  (4>) 

-sin  (0)  cos  (<J>) 
and  a  new  coordinate  system  by 
n  =  R  •  rT 


COS  (<f)) 

-sin  (0)  sin  (4>) 


sin  (0) 
0 

cos  (e) 


(16) 


(17) 


(18) 


(19) 


(20) 


where : 


4 


(21) 


The  new  coordinate  system  has  the  advantage  that  n  =  n  and  n  =  n 

yt  yo  zt  zo 

Only  nx  varies.  The  concentration  is  given  by 


C  = 


2QAQ 


(2tt)3/2o  ct  o 
x  y  z 


exp 


-n  n 


The  concentration-pathlength  is  then  given  by 


(22) 


sin2  (6)^  1 


2  ill  fXt 

exp  -l/2|ny+nz}  J  exp 


which  is  equivalent  to  Equation  12  if  the  square  of  the  argument  is  completed 
and  the  variable  of  integration  shifted.  The  value  of  Equation  23  is  that  it 
may  be  written  in  terms  of  error  functions  (erf's). 

The  error  function  is  given  by  ^ 


dn 


(23) 
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erf (x) 


(24) 


Equation  23  may  be  rewritten  in  this  form  by  defining  a  new  variable  of 
integration 


and  simple  algebra.  This  gives 


(25) 
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(26) 


where : 


sgn(x)  =  1,  x>0, 

=  1,  x-.O,  and  (27) 

=  0,  x=0. 


The  use  of  the  sgn  function  arises  from  the  fact  that  the  integrand  of  the 
error  function  is  even  (non-negative)  so  that  the  sign  of  the  integral 
depends  only  on  the  sign  of  the  end  point. 


Equation  26  is  usually  the  most  efficient  form  to  use  in  computing 
concent  rat ion-pathlength  since  the  error  function  may  be  compactly  computed 
using  a  single  Pade  approximant .  12  Alternately,  if  n  pxq  <  0,  |  nxt  |  >>  |, 

and  | Hxo |  >y\  >  then  the  integral  of  Equation  23  may  be  analytically  solved 
for  infinite  limits. H 


III.  APPROXIMATION  TO  THE  ERROR  FUNCTION 

A  function  that  is  well-known  to  solid  state  physicists  and  statistical 
mechanicians  is  that  used  to  describe  the  distribution  of  states  of  particles 
obeying  Fermi-Dirac  statist ics . ^  This  function  has  the  form 

f  (x)  =  [1  +  exp  (x)  ]  1  (2d) 

which  has  the  property: 
f (x)  =  1,  x  «  0, 

*  0,  x  »  0,  (29) 

=  1/2,  x  =  0. 

For  simplicity,  f(x)  is  referred  to  as  a  Fermi  function. 

Less  well-known  are  the  properties  of  the  derivative  of  the  Fermi  func¬ 
tion.  The  derivative  may  be  calculated  in  straightforward  manner  as 

df(xi  9 

—y - —  =  -exp  (x)  f (x)“  (30) 


=  -f(x)  f(-x). 
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The  derivative  of  the  Fermi  function  has  a  maximum  value  of  1/4  at  x  =  0.  It 
is  also  approximately  gaussian  in  shape.  This  is  shown  in  Figure  1. 

This  shape  leads  to  the  suspicion  that  the  Fermi  function  may  be  used  as 
an  approximation  to  the  error  function 


exp 


[V 


a  df (ax) 
d  (ax) 


(31) 


where  a  is  a  parameter  chosen  to  give  agreement  with  the  error  function. 


IV.  PERFORMANCE  OF  THE  APPROXIMATION 

As  an  investigation  of  the  relative  accuracy  of  the  approximation  for 
the  error  function  developed  in  the  last  section,  Equation  35,  the  SEMM 

model  ^  was  exercised  for  several  WP  munitions.  HC  (hexachloroe thana ) 
munitions  were  not  exercised  because  of  the  additional  computational  burdc  'l 
involved.  The  approximation  is,  of  course,  still  as  valid  for  HC,  and  re¬ 
duces  the  computational  burden,  but  the  time  dependence  of  the  calculation 
for  HC  only  detracts  from  the  examination  of  the  approximation. 

Several  geometries  of  target  and  observer  were  considered,  both  symmetric 
and  unsymmetric  with  respect  to  the  cloud.  The  quantity  of  comparison  used 
was  the  visual  transmission  through  the  cloud,  Equation  1,  between  observer 
and  target.  The  worse  cases  found  were  when  both  target  and  observer  were 
on  the  same  side  of  the  cloud  relative  to  the  cloud  centroid.  An  example  of 
these  results  are  shown  in  Figure  2  through  8. 

These  figures  are  plots  of  transmission  for  a  4.2-in.  WP  round,  neutral 
meteorology,  and  head  wind.  The  observer  and  target  are  located  on  the  x 
axis.  Thus,  the  LOS  either  passes  through  the  centroid  or  may  be  extended 
through  it.  The  observer  is  always  located  from  the  centroid.  The  dis¬ 

tance  between  observer  and  target  is  varied  from  0  to  6a^. 

It  may  be  seen  that  the  error  is  greatest  for  early  times  when  the  ob¬ 
server  and  target  are  within  of  each  other.  At  time  t  =  0,  this  distance 

apart  (a^)  is  about  4  meters,  increasing  approximately  linearly  with  time. 

The  error  decreases  with  time.  For  reasonable  observer  target  distances  of 

3  3  3 

10  m,  approximately  5  x  10  seconds  are  required  for  to  grow  to  10  m  in 

size.  This  is  interpreted  to  indicate  that  the  error  in  this  approximation 
is  not  stressing  for  most  simulation  uses. 


By  substituting  Equation  31  into  Equation  24,  we  may  get 


By  using  erf  (<»)  =  1,  Equation  33  may  be  written  as 


i  =  JL_  f(0) 


av  TT 


which  yields 


and 

erf  (x)  =  1  -  2f  /  —  \  (3 

Vr) 

Equations  26  and  35  may  be  combined  to  yield  an  approximate  concentrat ion- 
pathlength  of 


The  advantage  of  Equation  36  lies  in  its  simple  evaluation. 


V.  SUMMARY  AND  CONCLUSIONS 

A  closed  form  integrable  approximation  for  the  error  function  has  been 
developed.  This  approximation  has  the  advantage  that  it  can  significantly 
reduce  the  computat ional  burden  in  existing  obscuration  models  for  the  cal 
culation  of  concentrat ion-pat hlength.  This  computational  burden  currently 


represents  a  large  fraction  of  the  computer  memory  and  execution  time  required 
to  exercise  these  models. 

The  error  implicit  in  this  approximation  has  been  shown  to  be  such  that 
it  is  not  excessive  for  most  geometries  that  are  of  interest  to  obscuration 
modelers. 
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Figure  1,  Comparison  of  gaussian  and  derivative  of  Fermi  function. 
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Figure  Z  Comparison  of  transmission  calculated  using  error  function  and  Fermi 
function  through  smoke  cloud  at  t  =  0. 
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Figure  3.  Comparison  of  transmission  calculated  using  error  function  and  Fermi 
function  through  smoke  cloud  at  t  =  30. 
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Figure  6.  Comparison  of  transmission  calculated  using  error  function  and  Fermi 
function  through  smoke  at  t  =  240. 
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Figure  7,  Comparison  of  transmission  calculated  using  error  function  and  Fermi 
function  through  smoke  at  t  =  480. 
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Figure  8.  Comparison  of  transmission  calculated  using  error  function  and  Fermi 
function  smoke  at  t  -  960. 
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